FERMILAB-PUB- 13-082-E 



Measurement of the cross section for direct-photon production in association 
with a heavy quark in pp collisions at ^/s = 1.96 TeV 
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We report on a measurement of the cross section for direct-photon production in association with 
a heavy quark using the full data set of ^/s = 1.96 TeV proton-antiproton collisions corresponding 
to 9.1 fb~^ of integrated luminosity collected by the CDF II detector at the Fermilab Tevatron. The 
measurements are performed as a function of the photon transverse momentum, covering photon 
transverse momentum between 30 and 300 GeV, photon rapidities < 1.0, heavy-quark-jet 
transverse momentum p^** > 20 GeV, and jet rapidities ly^"*] < 1.5. The results are compared with 
several theoretical predictions. 
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collisions provides valuable information on the prob- 
ability distributions of partons inside the initial- 
state hadrons. At photon transverse energies 
[l[ smaller than 100 GeV, such events are produced 
predominantly by the Compton scattering process 
gQ jQ, while at higher energies the dominant 
process is quark-antiquark annihilation with gluon 
) splitting to heavy quarks qq ^ iQQ 

. It is conventional to assume that the charm 
(c) and bottom (b) quarks in the proton arise only 
from gluon splitting. However, there are other mod- 
els that allow the existence of intrinsic heavy quarks 
in the proton A cross-section measurement of 
'J + Q + X {X can be any final-state particle) pro- 
duction provides information on the heavy-quark 
and gluon parton distribution functions (PDFs) and 
on the rate of final-state gluon splitting to heavy 
quarks. 

The CoUider Detector at Fermilab (CDF) Col- 
laboration studied the process pp ^ j + b + X a.t ^/s 
= 1.96 TeV, for photons in the range 20 < < 70 
GeV [J]. The measured cross section agreed well 
with a prediction based on a perturbative quantum 
chromodynamics (pQCD) expansion Q at next-to- 
leading order (NLO) in the strong coupling constant 
as- The DO Collaboration measured the cross sec- 
tion of photons in association with heavy-flavor jets 
using data collected at -/s = 1.96 TeV, covering the 
range 30 < £;J < 300 GeV The results disagreed 
with the NLO pQCD prediction for both bottom jets 
and charm jets in the region > 70 GeV. 

In this Letter, we present the updated CDF 
measurements of the cross sections of photon with 
heavy-flavor jets, using the full data set from 9.1 
fb~^ of integrated luminosity collected by the CDF 
II detector, exploring up to 300 GeV with im- 
proved techniques. 

The CDF II detector @ has a cylindrical ge- 
ometry with approximate forward-backward and az- 
imuthal symmetry [ij. It contains a tracking sys- 
tem consisting of silicon microstrip detectors and 
a cylindrical open-cell drift chamber immersed in 
a 1.4 T magnetic field parallel to the beam axis. 
The silicon subsystem is used for reconstructing 
charged-particle trajectories (tracks) and heavy- 
flavor-decay vertices displaced from the primary in- 
teraction point. Pointing-tower- geometry electro- 
magnetic (EM) and hadronic calorimeters surround- 
ing the tracking system measure particle energies. 
At a depth approximately corresponding to the 
maximum development of the EM shower, the EM 
calorimeters contain fine-grained detectors (central 
electromagnetic strip chambers) that measure the 
shower profile. Drift chambers and scintillators lo- 



cated outside the calorimeters identify muons. 

The data are collected using a three-level online 
event-filtering system (trigger) that selects events 
with at least one energy cluster consistent with a 
photon in the final state. The trigger is approx- 
imately 100% efficient for signal events in the ex- 
plored kinematic region. 

The offline event selection requires the primary 
vertex z position to be within 60 cm of the center 
of the detector. Each event is required to have at 
least one photon candidate that has pseudorapid- 
ity [H in the fiducial region of the central calorime- 
ter (approximately |ry| < 1.04). The transverse en- 
ergy of the photon is corrected to account for non- 
uniformities in the calorimeter response, and cali- 
brated using electrons from reconstructed Z boson 
decays. Photon candidates are required to have 
E^ > 30 GeV and to satisfy pre-selection require- 
ments on calorimeter and tracking isolation and 
the ratio of the energy measured in the hadronic 
calorimeter to the EM energy, as described in Ref. 
Q. An artificial neural network (ANN) is con- 
structed from calorimeter- and strip-chambers-shape 
information and isolation variables to further reduce 
background Q . The photon candidates are required 
to pass a suitable threshold on the ANN output 
(0.75) for optimal signal-to-background discrimina- 
tion. 

At least one jet must be present in each event. 
Jets are reconstructed using the JETCLU algorithm 
with a cone radius R = \/iA(t))'^ + (A?])^ = 0.4 
in the azimuthal angle and pseudorapidity rj space 
[l|. We select jets that have Et > 20 GeV and \ri\ < 
1.5. At least one jet is required to be classified as a 
heavy-flavor jet using a secondary- vertex tagger [3]. 
This tagging algorithm exploits the long lifetime of 
hadrons containing 5 or c quarks, and is based on the 
reconstruction of a displaced, or secondary, vertex 
using the reconstructed tracks. If multiple-tagged 
jets are present in the explored kinematic region, the 
one with the highest Et is selected. The selected jet 
is required to be reconstructed in a volume outside 
an rj (f) cone of R = 0.4 surrounding the photon 
candidate. 

After all the selection requirements, 214 336 
events remain in the data sample. Two main back- 
ground sources contribute to these events: jets 
misidcntificd as photons (false photons) and light- 
flavor jets mimicking heavy-fiavor jets. To estimate 
the rate of false photons, the photon ANN distribu- 
tion in data is fitted to a linear combination of tem- 
plates for photons and jets, obtained from a simu- 
lated inclusive photon sample using SHERPA [i3| and 
a dijet sample using pythia pjj], after applying all 
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the photon selection criteria except the requirement 
on the ANN output. A fit is performed in each 
interval, yielding prompt photon fractions (purities) 
between 77% and 94% in the ANN signal range. The 
resulting photon purities and one example fit are 
shown in Fig. [TJ The systematic uncertainties on 
the photon purities are estimated by varying the in- 
put variables to the ANN within their uncertainties. 
The dominant uncertainty on the shape of the ANN 
originates from the modeling of calorimeter isolation 
energy. The overall uncertainty is estimated to de- 
crease from 6% at = 30 GeV to 2% for > 70 
GeV. 
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between 30 and 300 GeV. The systematic uncertain- 
ties range from 15% to 30% and are dominated by 
the uncertainties in the simulated MsecVtx template 
shapes originating from the uncertainty in the mod- 
eling of tracking-system efficiency. Figure [5] shows 
the result of the fit for Enj^ between 40 and 50 GeV, 
as an example. 
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FIG. 2: Distribution of the secondary-vertex mass of 
tagged jets after applying the full selection, for photon 
candidates with 40 < -EJ < 50 GeV. The points are 
data, and the stacked, shaded histograms represent the 
estimated contributions from the fit of the 6-, c-, and 
light-quark jets and false photon background. 



FIG. 1: Photon purity as a function of for events 
restricted to the ANN signal region. The fit to the ANN 
distribution for photon candidates passing pre-selection 
requirements and with between 40 and 50 GeV is 
shown in the inset. 



Backgrounds to heavy-flavor jets arise from 
light-flavor jets where random combinations of 
tracks mimic a displaced vertex. The fractions of 
h- and c-jets are determined by fitting the invariant 
mass (MsccVtx) of the system of charged particles, 
assumed to be pious, originating at the secondary 
vertex, using the templates for 6-, c-, and light-quark 



jets constructed with pythia The contribu- 

tion to the AfsccVtx distribution from events with a 
false photon is modeled using dijet data, where one 
jet is required to deposit most of its energy in the 
EM calorimeter to mimic a photon and the other jet 
is required to pass all the heavy-flavor-jet selection. 
The loose photon requirement selects predominantly 
false photons. This background component is then 
constrained to the number of false photons from the 
ANN fits. After subtracting the contribution from 
events with false photons, 22% to 37% of the ob- 
served tagged jets are 6-quark jets, and 16% to 24% 
of the observed tagged jets are c-quark jets for E"^ 



section as a func- 
as da^^"'^ jdET^ = 



The differential cross 
tion of Elp is defined 
iV/-y/h(f,)/(Ai?JetrigeuF'C), where N is the number 
of data events in a given bin after applying 
the full selection, /-^ is the photon purity, /^(c) 
is the 6-jet (c-jet) fraction in events with true 
photons, Aii^J is the i?J bin size, e^p^g is the 
trigger efficiency, and £ is the integrated luminosity. 
The bin-by-bin unfolding factor euF combines 
corrections for acceptance, efficiencies of the photon 
selection and tagging algorithm and resolution 
effects to infer the results at the hadron level, using 
prompt-photon events simulated with SHERPA 10]. 
The numerator of the unfolding factor is obtained 
by applying the same requirements to the SHERPA- 
simulated events as the ones applied to data. The 
denominator of the unfolding factor is obtained by 
applying the same kinematic and isolation selection 
on the generated quantities. Unfolding factors 
obtained with pythia [lH are used to evaluate the 
systematic uncertainties. The photon efficiency is 
calibrated by comparing the selection efficiencies for 
Z'^ — !■ e~ events in data and in simulation. The 
tagging efficiency is calibrated with data enriched 
with heavy-flavor jets. The unfolding factors range 
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from 18% to 27% for 7 + 6 + X events and from 
4% to 8% for 7 + c + X events. The systematic 
uncertainties are estimated to be approximately 
10% and are dominated by the uncertainties in the 
photon-energy scale and the tagging efficiency. 

The differential cross sections for 7 + 6 + X and 
') + c-\'X productions are extracted in eight intervals 
of i^J, and are listed in Table U and shown in Fig.[3l 
The total systematic uncertainties range from 18% 
and 35% and are dominated by the uncertainties in 
the 6-jet and c-jet fractions. A 6% constant uncer- 
tainty comes from the Tevatron integrated luminos- 
ity. 

The predictions based on NLO pQCD Q are 
given in Table |T] and shown in Fig. [31 These predic- 
tions include direct-photon production subprocesses 
and subprocesses where the photon is emitted from 
parton fragmentation, both at 0{aal). The calcu- 
lation utilizes CTEQ6.6M parton distribution func- 
tions [1^. The scale dependence is evaluated by 
varying the renormalization, /i^, factorization, ^y, 
and fragmentation, fip, scales, assumed to be the 
same, from the default value pip to and 2p'^. 

The predictions based on a /cy-factorization ap- 
proach [IJI are reported in Table |T] and Fig. |3l This 
calculation is based on 0{aa^) off-shell amplitudes 
of gluon-gluon fusion and quark-(anti)quark inter- 
action subprocesses, and the fcr-dependent (i.e., un- 
integrated) parton distributions, where kr denotes 
the transverse momentum of the parton. The non- 
vanishing transverse momentum of the colliding par- 
tons leads to a broadening of the photon transverse- 
momentum distribution. The scale dependence is 
evaluated in the same way as the NLO calculations. 

Both the NLO and fcr-factorization predictions 
are parton-level calculations without modeling of 
underlying-event energy. We correct those two pre- 
dictions for the non-perturbative effects of parton- 
to-hadron fragmentation and for underlying-event 
energy, by multiplying with a correction factor de- 
rived from a sample simulated with SHERPA. The 
correction factors are shown in Table HI 

The predictions of SHERPA [l3| are reported in 
Table U and Fig. |31 This calculation includes all the 
tree-level matrix-element diagrams with one photon 
and up to three jets, with at least one & jet or c jet 
in the explored kinematic region. This calculation 
features a parton-jet matching procedure in order to 
avoid an overlap between the phase-space descrip- 
tions given by the fixed-order matrix-element sub- 
processes and the showering and hadronization in 
the multi-jets simulation. 

The predictions of pythia [ii[ are reported in 
Table I] and Fig. |3l This calculation includes the 



2 — 7> 2 matrix-element subprocesses gb — >■ jb and 
qq — ^ 73 with g ^ bb and g ^ cc splittings in 
the parton shower. In the ratio plots, we multi- 
ply the pythia calculations by an empirical factor 
of 1.4 to improve the a gree ment of the normaliza- 
tion. Previous studies [12| showed that the con- 
tribution of gluon splitting to heavy flavor has to 
be approximately doubled over expectations from 
the leading-order PYTHIA generator to reproduce the 
data. Hence, we also show predictions that include 
a double gluon-splitting rate to heavy flavors. 

The NLO pQCD predictions agree with data 
at low Ep but fail to describe data for > 70 
GeV for the bottom-jet cross section. The same 
trend is observed in the charm-jet cross section even 
though the experimental uncertainty is larger. For 
large E^, the dominant production process yielding 
a photon and a heavy quark involves a final-state 
gluon splitting into a heavy-flavor pair. This pro- 
cess is present only at leading order in the NLO 
calculation. The SHERPA prediction allows up to 
three partons in the final state, through the inclu- 
sion of additional tree-level amplitudes. The addi- 
tional amplitudes also serve as a source of heavy- 
flavor pairs (through gluon splitting), which is im- 
portant for the high E^ range. The fcT-factorization 
and SHERPA predictions are in reasonable agreement 
with the measured cross sections. The pythia pre- 
dictions disagree with the data both in rate and in 
shape. Scaling the pythia prediction and doubling 
the rate for g ^ bb or g ^ cc leads to an improved 
agreement with the data. 

In conclusion, we measure the differential cross 
sections for inclusive production of a photon in as- 
sociation with a heavy flavor quark for E^ between 
30 and 300 GeV using the full CDF Run II data set 
and compare the results with four theoretical pre- 
dictions. The results indicate that an improved un- 
derstanding of gluon-splitting rates to heavy flavors 
is important for the NLO pQCD calculations and 
the PYTHIA generator to model data. These results 
can be used to improve the background modeling 
in the searches for new physics in channels involving 
the production of photons in association with heavy- 
flavor quarks and to test the models that contain 
intrinsic heavy quarks. 
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FIG. 3: The measured differential cross sections compared with theoretical predictions. The left panels show the 
absolute comparisons and the right panels show the ratios of the data over the theoretical predictions. The pythia 
predictions are scaled by 1.4 in the ratio distributions. The comparisons are shown for 'y + b + X (top) and 7 + c + X 
(bottom) processes. The shaded area around the data points indicates the total systematic uncertainty of the 
measurement. The scale uncertainties are shown for the NLO and the fcT-factorization predictions. 
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[1] A cylindrical coordinate system, (r, 0, z) is used 
with origin at the geometric center of the detector; 
r is the radius from the nominal beam line, tfi is 
the azimuthal angle, and +z points along the di- 
rection of the proton beam. The polar angle 9 with 
respect to the proton beam defines the pseudorapid- 
ity 77, which is given by rj = — ln[tan(^/2)]. Trans- 
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TABLE L The j + b + X and 'y + c + X cross sections in intervals of together with statistical and systematic 
uncertainties. Four theoretical predictions are shown. The scale uncertainties are shown for the NLO and the kr- 
factorization predictions. The SHERPA and pythia predictions have large scale uncertainties, which are not shown 
in the table. The last column shows non-perturbative corrections applied to the NLO and the fcT-factorization 
parton-level predictions. 



E'^ bins 


dcr/d.Bj 


NLO [2] 


kr fact. [14] 


SHERPA I 


10] 


PYTHIA [11] Corr. 


(GeV) 


(pb/GeV) 


(pb/GeV) 


(pb/GeV) 


(pb/GeV) 


(pb/GeV) 




-y + b + X 


30-35 


1.47±0.07±0.41 


2.09 ±0.10 


1.76 ±0.64 


1.84 




1.09 




0.937 


35-40 


(8.90±0.49±2.49)xl0"^ 


1.16 ±0.08 


1.05 ±0.34 


1.16 




7.38x10" 


-1 


0.936 


40-50 


(4.87±0.25±1.26)xl0"^ 


(5.18 ±0.54) X 10"^ 


(4.89 ± 1.67) X 10"^ 


6.04x10" 


-1 


3.44x10" 


-1 


0.915 


50-70 


(1.60±0.09±0.40)xl0~^ 


(1.53 ±0.22) X 10~^ 


(1.60 ±0.51) X 10"^ 


2.08x10" 


-1 


1.02x10" 


-1 


0.966 


70-90 


(5.17±0.51±1.41)xl0"^ 


(3.59 ± 0.70) X 10"^ 


(4.24 ± 1.21) X 10"^ 


5.83x10" 


-2 


2.94x10" 


-2 


0.954 


90-120 


(1.79±0.18±0.50)xl0"^ 


(9.45 ±2.35) X 10~^ 


(1.25 ±0.30) X 10-2 


1.79x10" 


-2 


8.22x10" 


-3 


0.920 


120-170 (4.49±0.81±1.58)xl0"^ 


(1.98 ±0.59) X 10"^ 


(3.13 ±0.51) X 10"^ 


4.19x10" 


-3 


1.94x10" 


-3 


0.907 


170-300 (6.39±2.26±2.04)xl0"'' 


(1.90 ±0.67) X 10"* 


(3.99 ±0.25) X 10"* 


4.30x10" 


-4 


2.37x10" 


-4 


0.913 


-y + c + X 


30-35 


(1.16±0.05±0.20)xl0 


(1.74 ±0.10) X 10 


(1.07 ±0.66) X 10 


1.25x10 


8.01 




1.28 


35-40 


6.33±0.33±1.08 


8.82 ±0.72 


6.22 ±2.77 


7.23 




4.39 




1.25 


40-50 


2.92±0.17±0.48 


3.67 ±0.36 


2.65 ±1.67 


3.43 




2.01 




1.21 
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(7.62±0.60±1.39)xl0"^ 


(8.54 ± 1.03) X 10"^ 


(7.26 ±3.02) X 10"^ 


9.79x10" 


-1 


5.12x10" 


-1 


1.16 


70-90 


(1.67±0.35±0.37)xl0~^ 


(1.62 ±0.25) X 10"^ 


(1.71 ±0.54) X 10"^ 


2.28x10" 


-1 


1.05x10" 


-1 


1.13 


90-120 


(4.37±1.44±0.85)xl0"2 


(3.51 ±0.65) X 10-2 


(4.99 ± 0.97) X 10"2 


5.90x10" 


-2 


2.50x10" 


-2 


1.11 


120-170 (1.32±0.55±0.26)xl0"2 


(5.44 ± 1.37) X 10"^ 


(1.25 ±0.02) X 10-2 


1.20x10" 


-2 


4.56x10" 


-3 


1.07 


170-300 (1.51±1.23±0.45)xl0"^ 


(3.86 ± 1.16) X 10"'' 


(1.92 ±0.10) X 10"^ 


1.12x10" 


-3 


4.84x10" 


-4 


1.04 
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